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Abstract 

We obtain an identity between Fredholm determinants of two kinds of operators, one acting 
on functions on the unit circle and the other acting on functions on a subset of the integers. 
This identity is a generalization of an identity between a Toeplitz determinant and a Fredholm 
determinant that has appeared in the random permutation context. Using this identity, we prove, in 
particular, convergence of moments for arbitrary rows of a random Young diagram under Plancherel 
measure. 

1 Introduction 

In jjj, the authors considered the length In{k) of the longest increasing subsequence of a random 
permutation n G Sn, the symmetric group on N numbers. They showed, in particular, that for 
«Jv(tt) := — ^^tts , 

lim F(£ N < x) = F {1) (x), (1.1) 

where F^\x) is the Tracy- Widom distribution [|3(| for the largest eigenvalue of a random matrix from 
the Gaussian Unitary Ensemble (GUE). The authors also proved the convergence of moments, 

/CO 
x m dF^(x), m = l,2,--.. (1.2) 



The authors then reinterpreted (1.1), (1.2) in terms of Young diagrams A = (Ai,A2,---) via the 
Robinson-Schensted correspondence. Here Aj is the number of boxes in the jth row of A and Ai > 
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\ 2 > . . . > o. The set of Young diagrams Yjv of size N, y\ Xj = N, is equipped with Plancherel 
measure, 

P^ lan (A) := , A e Y N , (1.3) 
where d\ is the number of standard Young tableaux of shape A. Set 



Then (LI), (L2) imply that £i converges in distribution, together with all its moments, to This 
reinterpretation led the authors to conjecture that for all k, £i, £21 " ' 1 £fc converge to the joint distribu- 
tion function F(x\, X2, • • • , x k ) for the first k eigenvalues of a random GUE matrix. In the authors 
verified the convergence in distribution, together with its moments, to the Tracy- Widom distribution 
for the second largest eigenvalue of a random GUE matrix. The conjecture for £1,^2, ■ • • ,6e was 
then proved in three independent papers [^9| , || , |24| , all appearing within a few months in the spring 
of 1999. Let yj be the jth largest eigenvalue of a random N x N matrix from GUE with probability 
density 

1 N 
dF%™( yi ,... ,y N ) = — [] (y l -y ] ) 2 H^dy 1 ---dy N , (1.5) 

N l<i<j<N j=l 

where y± > • • • > yjy, and Zjq is the normalization constant. At the 'edge' of the spectrum, the following 
convergence in distribution is well-known (see, e.g. [M Theorem 1.4): for any k € N, there is a 

distribution function F(x±, • • • , CEfc) on 11 > • • • > Xk such that 



lim P^ UE (( yi - V^V)^ 1 / 6 <x lr -- ^yk-Vw^N 1 / 6 < x k ) =F(x ir - ,x k ). (1.6) 
In all three papers |29[, S, p4|, the authors showed that for any xi, • • • ,Xk € 



hm P£ lon (£i <x lr -- <x fc ) = F(xi,--- ,x fe ), (1.7) 

TV— > 00 

but the question of the convergence of moments was left open. 
Introduce the Poissonized Plancherel measure 

Pf- (A )=£^_P«««(A), t>0, (1.8) 

AT=0 

on all Young diagrams, which corresponds to choosing TV as a Poisson variable with parameter t . 
Here P^ lan (A) = if A is not a partition of N. Throughout the paper, we will work with Pf ois (A) 
rather than P^/ an (A) itself. This is because the expectation with respect to Pf OIS (A) leads to convenient 
determinantal formulae. Indeed, in |lS|| , Gessel proved the following formula 

Pf ols (Ai < n) = e- e det(T n ), (1.9) 
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where T n is the n x n Toeplitz matrix with entries (T n ) P9 = c p _ 9 , < p, q < n, where c k is the k 



Hi 



Fourier coefficient of e*( z+z 



J\z\=i- 



-k e t(z+z )-^|-, This formula played a basic role in || in 



proving (1.1), (1.2). In Q, the authors introduced the integral operator K n with ip(z) = e tlyZ z 
( |2.l|) below) and proved the following formulae 



$Pois 



(Ai <n) = 2" n det(l - K n ) 



(see 



(1.10) 



and 



(A 2 <n+l) 



(Ai < n) + 



_5 



[(l + ^-"det(l-^K n )]. 



(1.11) 



These formulae played a basic role in || in proving the analogue of (1.1), (1.2) for A 2 . In [g| and |24|] , 
and also later, in greater generality, in |28| and fl31| , the authors obtained the following identity : Let 
Afe denote the (finite) set {n £ {0, 1, • • • } k : Ysj=i n i — r ~ 1' r = ■"•> ' ' ' > Then for a& < • • • < a% < 

ao = oo, 



Pf ois (Ai 



nGA) 



1 < ai, A 2 
1 



2 < a 2 , 



, Afc - k < a k ) 



■ n k \ ds" 1 



■■ ds T 



k 
\ " 



det(l + 

si=—=s fc =-l i=i 



(1.12) 



where the matrix elements of S(i, j) are given in ( pT3| ) below with (p(z) — e^( z ~ z ± \ As usual, X(a,b] 
denotes the characteristic function of the interval (a, b], and so (X^=i s iX(ai.ai_i]) S denotes the operator 
in ( 2 (Z) with kernel s; S(i, j) if i € (a;,a/_i], and zero otherwise. Setting = It + Xjt 1 ^ 3 , xi > x 2 > 



> Xfc, and letting t — ► oo, and de-Poissonizing as in f2q| , the authors in [|8| and |24|] obtain (1/7). In 



and 



however, the authors are not able to prove convergence of moments. The reason for this is 



that it is possible to use the classical steepest-descent method to control det(l + (^=i s iX(a 



S) 



for o, = 2t 



Xjt 1 / 3 as t 



oo, uniformly for x\ > x 2 > • • • > x k > M for any fixed M. But as the 



Xj's tends to -co, the method break down. On the other hand, the authors in || |] are able to control 
the lower tails of the probability distributions, and hence prove the convergence of moments for Ai and 
A 2 , using the steepest-descent method for the Riemann-Hilbert problem (RHP) naturally associated 
with T n and K n above. The steepest-descent method for RHP was introduced in |Q, and extended to 
include fully non-linear oscillations in |Lf|. The asymptotic analysis in ||, Q is closely related to the 
analysis in |l^, |l4| . The main motivation for this paper was to find a formula for the joint distribution of 
Ai, • • • , Afc, which generalized (1.11), and to which the above Riemann-Hilbert steepest-descent methods 



could be applied to obtain the lower tail estimates. 



Note that from (1.9), ( 1. 1C ) and (1.12), we have three formulae for the distribution of Ai, 



?f oU {Xi <n) = e-* det(T n ) 

= 2 -B det(l - K n ) 
= det(l - X[n,oo) S) 



(1.13) 
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and from (1.11) and (1.12), two formulae for the distribution of A2, 



" ols (A 2 < n + 1) = Pf ois (A! <n)+[~ 



[(1 + Vi)""det(l- v^K n )] 

L 

det(l + sx[ n ,oo) S). 



(1.14) 



To obtain the second formula, we use the fact that Ak=2 — {(0, 0), (0, 1)} and set a\ — 00, = n — 1 



in (1.12). From (1.14), we might guess that 



(1 + yfs)- n det(l - K n ) = det(l - s X[n ,oo) S). 



(1.15) 



The content of Theorem 2.1 is that precisely this relation is true for a general class of functions f>[z), 
provided f(z) has no winding. If the winding number of ip is non-zero, the above relation must be 
modified slightly as in ( |2.7| ). The fact that e~* det(T n ) = det(l — X[n,oo) S) for (essentially) the same 
general class of tp's (with zero winding number) was first proved in Q, with an alternative proof given 
in |||. The relation (1.15) for general s was proved essentially simultaneously with the present paper by 
Rains in plj , for a subclass of functions (p with zero winding, using algebraic methods (see Remark 4 
in Section |^). A particularly simple proof of the relation e~* det(T n ) = det(l — X[n,oo) S) can be found 
in the recent paper g of Bottcher (see also Q). The paper M also extends Theorem |2.1| and |2 . 1 2| 



to 



the matrix case (see Remark 2.3 and 2. IS below). 

In this paper, we will prove a general identity between determinants of operators of two types : the 
operators of the first type act on functions on the unit circle, and the operator of the second type act on 
functions on a subset of the integers. Specializations of this identity have, in particular, the following 
consequences : 



(51) A proof of the convergence of moments for £1, • • • , (see Theorem 3.1) 

(52) An interpretation of F(xx, • • • , in ( P~7t| ) as a "multi-Painleve" function (see Section ^J). As we 
will see, the behavior of multi-Painleve functions has similarities to the interactions of solitons in 
the classical theory of the Korteweg de Vries equation. 



(53) The analogue of Theorem 3.1 for signed permutations and so-called colored permutations (see 
Section ^) 

(54) New formulae for random word problems, certain 2-dimensional growth models, and also the 
so-called "digital boiling" model (see Section [jj) 



The new identity is given in Theorem 2.1 in two closely related forms (2.7), (2.S). In (S1)-(S4), we only 



use (p?7D. 

As we will see, some simple estimates together with a Riemann-Hilbert analysis of det(l — ^/sKn) 



is enough to control the lower tail estimation of P t ols (A). The relation (1.15) generalizes to the multi- 
interval described in Theorem |2.12| in Section ||. 
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In Section 0, we prove the main identity (2.7), (2.S) in the single interval case, and also the identity 



( 2.54 ) in the multi- interval case. In Section |3j, we use (2/7) to prove the convergence of moments for 



random Young tableaux (Theorem 3.1). A stronger version of this result is given in ( p.2| ). Section |4| 
contains certain tail estimates, needed in Section ||. Various estimates needed in Section || for a ratio 
of determinants are derived in Section |s| using the steepest-descent method for RHP's. In Section ||, 
we introduce the notion of a multi-Painleve solution, and in Section ^ we prove various formulae for 
colored permutations and also discuss certain random growth models from the perspective of Theorem 



2.1 
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2 Fredholm determinant identity 

Let (p(z) be a continuous, complex-valued, non-zero function on the unit circle X = {z G C : \z\ = 1}. 
Define K n to be the integral operator acting on L 2 (T,,dw) with kernel 

K4z, ro ) ^" Z yf^ W - , (K n f)(z)= f K n (z,w)f(w)dw. (2.1) 

21TI(Z-W) J\w\ = l 

For a function / on E, its Fourier coefficients are denoted by fj, so that 

/(*)=£/^. ( 2 - 2 ) 

jez 

Let S be the matrix with entries 

= 5Z(v _1 )i+fcV-j-fc 5 h0 e z > ( 2 - 3 ) 

k>l 

and let R be the matrix with entries 

= '^2( l P~ 1 )i+k(P-j-k, U'GZ, (2.4) 

fc<0 

Let S n denote the operator X[n,oo) S acting on £ 2 ({n, n + 1, • • • }), 

(Sn/)(i) = £S a (i,i)/(j), *>n, (2.5) 

and let R n denote the operator Xf-oo.n-i] ^ acting on £ 2 ({- ■ ■ , n — 2, n — 1}), 

(R* /)(*)= E R *(h3)f(j), i<n-l. (2.6) 

j<n—l 
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Theorem 2.1. Let <p(z) be a non-zero junction on the unit circle satisfying Xjez \ < 00 > which 
has winding number equal to -For s € C otic? n G Z, K n , S n and R n are trace cZass on L 2 (T,, dw), 

£ 2 ({n, n + 1, • • • }) and ^ 2 ({' • • , n . — 1}) respectively, and we have 

det(l-sK n ) = (l + s) n+#(v) det(l-s 2 S n ), (2.7) 
= (l-s)-"- #(v) det(l-s 2 R n ), s^l. (2.8) 

Remark 2.2. Standard Banach algebra estimates show that if the winding number of ip is zero and 
X IjVjI < £/ien || logt^ljoo + (X Ull(l°8^)j| 2 ) < 00 • ?7ms is enough to prove that the first and the 
third terms in ( 1.13Q are equal for all such ip's (see j^J). In particular, by (1.13), (2.7) is true for all tp 



without winding and satisfying X \j<Pj\ < 00 > when 8 = 1. 

Remark 2.3. As noted by Bottcher ||/, Theorem 2.1 remains true in the case where tp{z) is an invertible 



N x N matrix, provided the exponent n + #(ip) is replaced in (2.7), (2.8) by Nn + #(det ip). The proof 
in the scalar case extends to N x N matrices, and we give no further details : the proof in Ji^J is different 
and uses Wiener- H op f factorization directly. 



For the proof of Theorem 2.1, we use the following basic properties of the determinant (see, e.g., 



[p3[). If A is a trace class operator on a Hilbert space H, \\A\\i = (trA* A) 1 / 2 denotes the trace norm. 

Lemma 2.4. (i). If A n is a trace class operator for each n and A n — > A in trace norm, then A is a 
trace class operator and det(l + A n ) — > det(l + A) as n — + 00. 

(ii). If A is a trace class operator, and B n and C n are bounded operators such that {B n )* and C n 
converge strongly to B* and C respectively, then det(l + C n AB n ) — > det(l + CAB) as n — > 00. 

(Hi). If AB and BA are trace class operators, then det(l + AB) = det(l + BA). 

(iv). Suppose C acts on £ 2 (Z) and has matrix elements (cij)i,jei- tfY^ijez, \ c ij\ < °°> then C is trace 
class and \\C\\i < J2i,jcz \ c ij\- 

Proof of Theorem \2. 4 Define the projection operators on the circle 

(P n f)(z) = J2fjZ j , neZ, (2.9) 

j>n 

and 

(Q n f)(z)= ]T /M ™ >0 > 
0<i<n 

(2.10) 

(Qnf)(z) = - ]T h Z '> n<0 ' 
n<j<Q 

with {Qof)(z) = 0. Thus in particular, we have P n = Pq — Q n . Let M g denote the multiplication 
operator 

(M g f)(z)=g(z)f(z). (2.11) 
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Direct calculation shows that 



K n = -Pa + M tp P n M v -i = (1 - P ) - M V {1 - P n )M v -i, 

S n =P n M v - 1 (l-P )M ip P n , 

R n = (1 - P n )M v -iP M v (l - P n ). 



(2.12) 



First, we show that K n , S n and R n are trace class. Indeed K n = —Q n — HM v -i, where H = [P n , M v \. 
H acts on the basis for £ 2 (£, dw), as follows : Hz k = J^i HikZ 1 . We find 

Hlk = <fl-k(Xl>n - Xk>n), l,k € Z, (2-13) 

where X>n denotes the characteristic function of the set {k > n}. But J^z k \^ik\ < J2j bVj'l < °°' an d 
hence by Lemma 2.4 (iv), we have the trace norm estimate 

II 111 <n + (53lM|W -1 ||£-. 



(2.14) 



Now write S n = AB where A : £ 2 ({1,2, ■ ■ ■}) -> l 2 {{n,n + l,---}) and B : ^ 2 ({n,n + l,---}) 
£ 2 ({1, 2, • • • }) with matrix elements 



^ifc = (<p 1 )»+*;, i > n, k > 1, 
-Bfej r = f-k-j, k>lj>n. 



Write 



^ = xt^X-iR 
where (Rf)j — f-j, Q" 1 denotes convolution on £ 2 (Z) by {(pj 1 }, 



(2.15) 



(2.16) 



(2.17) 



and Xn' X—i are the projections onto {k > n} and {fc < —1} respectively. From ( [2.16 ), it is clear that 
A is bounded from ^ 2 ({1, 2, • • • }) — > ^ 2 ({n, n + 1, • • • }) with norm estimate 



1141 < II^IU* 



(2.18) 



On the other hand, a similar calculation to (2.14) shows that B is trace class from £ 2 ({n, n + 1, ■ ■ ■ }) 
^({1,2,---}) and 



\B\\i< J2 i^-ii <El^' 



Z>n+1 



which implies 



s n ||i< ( £ K' + Mte- 



Z>n+1 



(2.19) 



(2.20) 
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Similarly, we have 



Rn||i< + \n\)v-i\) | 



Kn 



(2.21) 



Thus K n , S n and R n are trace class. Moreover, if we set tpj := Yl\j\<j L P3 Z ' > i J — 0, then from (the 
proofs of) Q2.14D , Q2.20Q and fl2.21| ), it is clear that as J -> oo, K n (<^j) -> K n (^), S n (v3j) S„(y>), 
R n ((^j) — » R n (ip) in trace norm, and hence the Fredholm determinants converge to the corresponding 
determinants. Also for J sufficiently large, the winding number of ipj is the same as the winding number 
of ip, and so we see that to prove (2.7), it is enough to consider <^'s which are non-zero and analytic in 
a neighborhood of S. Henceforth we will assume that ip is analytic : this analyticity assumption is not 
necessary and is used only to give a particularly simple proof of Lemma 2J3 below. 

In the below, we only present the proof of fl2.7p . The proof of fl2.8| ) is similar. Formally, we proceed 
as follows. Suppose P n is finite rank so that Pq = Q n + P n is also finite rank. We have 



det(l + sP - sM v P n M -i) = det(l + sP ) det(l - 



1 + sPo 



M v P n M v -,) 



Using P = Q n + P n and ,,\ p = 1 — tt-Po, the right-hand-side reduces to 



det(l + sQ n ) det(l + sP n ) det(l - s(l 



l+s 



P )M v P n M v -i 



n > 0, 



(2.22) 



(2.23) 



(det(l - sQ n )) 1 det(l + sP n ) det(l - s(l - -^- s P )M v P n M^i), n < 0. 

and P„ = P^ for the last 



The first term in both cases is equal to (1 + s) n . Using Lemma 2.4 
determinant, (2.23) becomes 

(1 + s) n det(l + sP n ) det(l - aP n M v -i{\ - -^—P )M v P n ) 



= (1 + s) n dct((l + sP n ) - a(l + s)P n M v -i (1 - -—P )M v P n ) 

l + s 



(2.24) 



= (l + s)"det(l-s 2 S n ), 

which is the desired result, up to the winding number #(y). For the case in hand, however, P is not a 
trace class operator and the above "proof" breaks down. We circumvent the difficulty by approximating 
the operator K n by finite rank operators, and the missing factor will appear along the way. 

Let Tjv be the projection 



(T N f)(z) = fi 



N > 1. 



|j|<JV 



Note that Tn is a trace class operator since it has finite rank. Clearly T^,T l 



N 



(2.25) 

1 strongly, and hence 



by Lemma 2A (rf) 



det(l - sK n ) = lim det(l + s(P - M v P n M -i)T N ). 



N 



(2.26) 



Now since PfeTW is trace class, proceeding as above in (2.22)-(2.24), we have for N > 



n , 



det(l + s(P - M^PnM^Tx) = (1 + s) n det((l + sP n T N ){l - sP n M v -iT N {l 



-P T N )M v P n )). 

i 

(2.27) 
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Thus we have 



det(l-sK n ) = (1 + s)" lim det(l + X N + Y N ), (2.28) 

N— >oo 



where for TV > 



X N = a 1 ** N P n (T N - M v -iT N M v )P n , (2.29) 

J. ~j~ s 

Y N = -s 2 1 + f N P n M v -iT N (l - P )M v P n . (2.30) 



We observe that 

(i) . Xn and Y/v are trace class. 

(ii) . X^ — > strongly as N — > oo. 

(hi). Yat — > — s 2 S n in trace norm as A^ — > oo. 

(iv). (1 + Xjv) -1 and (1 + Y/v)^ 1 are uniformly bounded in operator norm as N — > oo when s is small 
enough. 



The third property follows using Lemma 2.4 (ii) as (1 — Po)M v P n is trace class and T/v — ► 1 strongly. 

For a moment, we assume that s is small so that (iv) is satisfied. Now we rewrite the right hand 
side of fl2.28| ) as 

det(l + X N +Y N ) =det(l + y JV )det(l + X jV )det(l- - 1 1 Y N X N ). (2.31) 

1 + Ajv 1 + x JV 

From the properties (i), (ii), (iv) above, we have 1 , Xn 1+y Yn^n — > in trace norm. Using the 
property (iii), we now have 



det(l - s K n ) = (1 + s) n det(l - s 2 S n ) lim det(l + (2.32) 

N— »oo 



Rewrite Xjv as 



X n = s 1 ^f rN P n (M v -i(l - T N )M V - (1 - TAr))P„ 

= 8 ^* N P n (M 9 -iP N +lM v - P N +l)Pn + 8 ~^* N P„(M 9 -l(l - P- N )M V ~ (1 ~ ^-Jv))P„ 

=: Zjv + W N . 

(2.33) 

Then as AT > Inl, 



Wat = a 1 , . - P n M v -i(l - P- N )M v P n = s > „M r i(l - P_ N )(l - P Q )M v P n , (2.34) 



1 + sTm , , 1 + sTm 

—JLp n M v -i{l - P- N )M v P n = s—-^ 
1 + s * 1+s 

and hence Wm — * in trace norm as (1 — P-jv) - * strongly and (1 — Po)M v> P n is in trace class. Also 

Zn — > strongly, and (1 + Zn)^ 1 , (1 + VKjv) -1 are uniformly bounded as A~ — » oo for s small enough. 

Thus by similar arguments leading to fl2.32j ), we have 

lim dct(l + XAr)= lim det(l + Zjv). (2.35) 

AT— >oo N— >oo 
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Now by fl2.12| ) and wc note that 

M v -iPn + iM v — Pn+i = Kn+i +Qjv+i = An+i Kq A n+1 , 



(2.36) 



where Kn+i is Kn+i with tp replaced by (p 1 , and A/v+i is the operator of multiplication by z 
Thus we have 

det(l + Z N ) = det(l + s + S N P n A N+1 K A^ +1 P n ) 



JV+l 



1 + s 



= det(l + s K A N x +1 P n 1 ± * Tn P n A N+1 ), 

_L — \~ s 

by Lemma (iii). Since PkA^+i — Apj + iPk-N-i, and Tjv = P-n — P/v+i, we have 



(2.37) 



1 



\iV+l J 



1 + S 



1 + S 



((1 + S)P„_ 2V _ 1 - sP ) 



1 + s - sPp 
1 + s 



D(s) 



strongly. Since Kq is a trace class, we obtain 



lim det(l + Z N ) = det(l + sK Q D(s)). 



Therefore from (|3|) and fl2.35| ), 

det(l-sK n ) = (l + s)™dct(l-s 2 S*„)det(l + sK L>(s)). 



(2.38) 



(2.39) 



(2.40) 



Since det(l + sKq D(s)) does not depend on n, we obtain the value of this determinant by letting 



n — > oo in both sides of (2.40). But by Lemma [D] below, for small s, (1 + s) "det(l — sK n ) — > 
(1 + s)*^ as n — > oo. On the other hand, from ( 2.2C ), det(l — s 2 S n ) converges to 1. Therefore 



det(l + sK D(s)) = (l + s)*&\ and we obtain, for small s, 

det(l - s K n ) = (1 + S ) n+ * M det(l - s 2 ^), 
as desired. The result for all s now follows by analytic continuation. 



(2.41) 
□ 



Observe from (|+]), j2.12j ) that || K n || < 1 + | (^||oo||y 1 ||oo —'■ so- (A different estimate (see Appendix 



) shows that ||i^n|| < max (IMI 



} ).) Then we have the following result. 



Lemma 2.5. For a function ip which is analytic and non-zero in a neighborhood of the unit circle 
{\z\ — 1} in the complex plane, and has winding number equal to #(</?), we have for \s\ < Sq , 



lim (1 + s)~ n det(l - s K„) = (1 + s 



,#(¥>) 



(2.42) 



lim (1 - s) n det(l - sK n ) = (1 - s)-* {lp \ 



n - x 



(2.43) 



for some sq > 0. 
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Proof. In Lemma 5 of we obtained the result ( |2.42| ) for p(z) = e ^( z - z ) and s = < t < 1. For 
general analytic p and s small, the proof remains the same until equation (50). The second component 
in the asymptotics of F' is now — jfr^r, and hence (51) is changed to 



logdet(l - sK„) 



r ds' 




J 2m(l + s') 


j 

J\z\ = ll 



-nz 1 — 



p'{z) 



dz + 0(e- cn ) 



(n + #(tp))\og(l + s) + 0(e- cn ). 



The calculation for (2.43) is similar and we skip the details. 



(2.44) 



□ 



Remark 2.6. Lemma 2.5 does not require p to be analytic. However, in this case, the proof is particularly 
simple, and can be quoted directly from W as above. 

Remark 2.7. The fact that det(l + sK D(s)) = (1 + s)*^ is rather remarkable. It is an instructive 
exercise to check this identity directly when ip is simple, say p = z k or p — (l + az)(l + bz~ 1 ) , \a\, \b\ < 1. 



Remark 2.8. By (2.7), we see that ifn + #(ip) > 0, det(l — sK n ) has a root at s = —1 of order at least 



n + f/:(p). In particular, K n has eigenvalue —1. Moreover, i/K n is self-adjoint (which is true by (2.12) 
whenever \p\ — 1, e.g., p — e ^( z ^ z ) as in j^, Q/j, then K n has an eigenspace of dimension at least 
n + ffc((p) corresponding to the eigenvalue — 1. It is also clear from ( |2.7| ) that if s ^ —1 is a root of 
dct(l — s K n ), then so is —s. On the other hand, if n + < 0, then clearly det(l — s 2 S n ) has a root 
at s — ±1, etc. In the self-adjoint case, when \p\ = 1, we see from ( |2.12 ) that S n is positive definite 
with norm < 1 . We will use this fact in Section [|. 



Remark 2.9. Define the operator A acting on £(!*) by 



(2.45) 



(2.46) 
(2.47) 



Since det(l — sK n ) = det(l — sM v -i K n M v ), using ( [2.12 ), the above theorem can be rephrased as 

det(l - s(P n - A)) = (l + s) n+ #^det(l-s 2 P n (l-A)P n ), 

= (1 - s )- n -#^ det(l - s 2 {\ - P n )A{\ - P n )). 

These are the identities (8.55), (8.56) in for a certain subclass ofp's with zero winding, #(<p) = 

The following Corollary will be used in the analysis of (S3) in Section [t] below. 
Corollary 2.10. Let p(z) be as in Theorem 2.1. Define S^ 111 -* and R( m ) to be the operator; 



analogous to K n , S*" 1 ) and R/ m ) with the matrix elements given by 

1 - z n p{z m )w- n p(w m Y 1 



2m(z — w) 



(2.48) 
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and 



S (m) («,j) = 1 )( l +/c)/m^(-i-fc)/m, 



(2.49) 
(2.50) 



fc<0 



where tp a — (ip 1 ) a = i/ a ^ Z. S'ei si'"'' = X[n,oo) S^ m ^ and Ri" 1 ' = X(-oo,n-i] R^ m '. 27ien we have 



(m) 



det(l- sKW) = (l + s)" +# ^)det(l-s 2 S( 1 m) ), s ^ -1, 
= (1 - fl )-"-#(*0 det(l - s 2 4 m )), s ^ 1. 



Remark 2.11. Observe that s["^ has the block structure 



( Sn" 1 ^ (mi, mj) 

(m) 



Si m ^(mi, mj ' + to — 1) \ 



ySn' 1 ' (mi + TO — 1, mj) • ■ ■ Sn ( m * + TO — 1, TO J + TO — l)y 



/Sn(*,i) o 

S„(i,j) 
\ 



(2.51) 
(2.52) 



\ 




Sn(i,i)/ 

(2.53) 



For the multi-interval case, we can generalize the argument in Theorem 2.1 to obtain the following 
result. 

Theorem 2.12. Let Q = n$ < n\ < n?, < ■ ■ ■ < flu < nk+i — oo be integers, and let Si, • • • , sjt &e 
complex numbers satisfying Sk 7^—1 and Sfc — Sj ^ — 1. ylfeo set sq = 0. WTe ftawe 



det(l-^(sj -Sj_i)K nj ) 

rfe-l 



(2.54) 



= (1 + Sfe )#(*>) 
where #(tp) is again the winding number of ip 



H(i + s k - Sj ) n ^-^ 

3=0 



det 1- El 



"X[nj,n j+ i) I S 



Remark 2.13. As noted by Bottcher (Mj ; ef. Remark 2.S above), the formula (2.54) remains true in 
the N x N matrix case, provided we replace — Tij by N(nj + i — rij), < j < k — 1, and #((p) by 
#(det<p). Again the proof in the scalar case extends to the matrix case, and we provide no details. 

Proof. The formal procedure (without considering the winding number) is as follows. For j = 0, • • • , k — 
1, let Rj be the projection operator on {nj, ■ ■ ■ , %+i — 1}, and let Rk be the projection operator on 
{nk, nk + 1 • ■ • }• Since we have from ( [2.12j ) 



K„ =-^Ri + M tp [Y,Rl)M v -i, j = l,---,k, 

"1=3 



(2.55) 
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the determinant on the left-hand-side in (2.54), denoted by (*), is equal to 



(*) = dct ( 1 + s k Rj -MJY, SjRj ) m v 

3=0 S=l 



(2.56) 



First we pull out the term 1 + s k £ j= q R i' then use Lemma [2.4| (hi) to obtain 



(*) = dct ( 1 + s k ^2 R j ) det ( 1 



3=0 
k 



3=0 



= det I 1 + s fe ^2 R 3 ) det ( 1 _ ( X s * R i ) M ' 



3 = 1 



1 



(2.57) 



M 

1 + Sk £ i=0 R 3 



Now note that (recall sq = 0) 



3=0 



det ( 1 + s k y^Rj ) = dct I 1 + y^(s fe - Sj 

fc-i 



fe-i 



3 = 1 J 

det ( 1 + ^](sfc - s^O-Rj ) det ( 1 



^SjRj) (2.58) 



A— 1 



J](l + Sfe - Sj .)»^-% detfl + ^ T 



3=0 



-R 



Using (2.58) and then multiplying two determinants, we have 



(*) = 



fc-i 



3=0 



det i 1 + £ - + *) (t IT^) 



Finally, using 



1 + Sfe 



in the determinant on the right-hand-side of (2.59), we obtain 



(*) 



fc-i 



3=0 



det ( J - E 1+ T k j - Sj R ^- 1 [ l -£ R ) M -) 



(2.59) 



(2.60) 



(2.61) 



which is precisely ( 2.54 ) from ( [2.12| ). 



The rigorous proof is also similar to the proof of Theorem 2.1. Let Tjv be the projection on \j\ < N 



as in ( 2.25 ). We take N large so that N > n k . The analogue of ( 2.28| ) is now 

rfe-l 



(*) 



H(l + s k - Sj ) n ^-^ 



3=0 



lim det(l + X N + Y N ), 



(2.62) 
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where 

u 

sj 1 + s k - Sj(l - T N ) 

-f 1 + Sk - Sj 1 + s k 



X 



N 



fe -, /-, m n \ k 

Y N = -(' 

\ j=1 - . ~ K ~j - , ~« / j=0 



E i .?/ J ' /V) ^)^, /v(l - £ S ^)M„ (2.64) 



which becomes, by the same argument leading to ( 2.35 ), 

rfc-l 



(*) = 



■j=0 



k 

Sk Sj 



+ s k ~ Sj r^- n * det 1 - £ fc 3 X K , % . +1) S Km det(l + Z N ), (2.65) 
v j~ i ' v 



with Z^v in (2.33) where s is replaced by s k , This then leads to the desired result as in the single interval 



case. □ 

3 Convergence of moments 

In this section, we prove the convergence of moments for arbitrary (scaled) rows, £j , of a random young 
diagram under the Plancherel measure, mentioned in the Introduction. The tail estimates used in the 
proof of Theorem [O] are given in Section [| below. 
Let N := NU {0}. 

Theorem 3.1. For any fixed k € N, and for any a,j € No, 1 < j ' < k, we have as N — ► oo 

^(z? ■■■%") ^Hxr ■■■x a k k ) (3.i) 

where E^ lan denotes the expectation with respect to the Plancherel measure on Yjy, and E denotes the 



expectation with respect to the limiting distribution function F in (1.6), (1.7) 



Remark 3.2. It will be clear from the proof below that the following stronger convergence result is also 

true : Let hj(x), j = 1, • • • , k be continuous functions on K satisfying \hj(x)\ < C\e C2 ^ 1 for some 
e > 0. Then for any k, as N — > oo, 

E^ lan (M£i) ■ • ■ -» • ■ • h k (x k )). (3.2) 

Proof. We have 

jgPlan^ ...£*)=/" 7J ^'dP^l < Si, • • • , Cfe < X k ) (3.3) 

Jxi> - >Xk j_i 

since Ai > A2 • • • • Fix a number T > 2. We split the integral into two pieces : 

(a) max IxA < T (3.4) 

i<j<fe J 

(b) max |x,-| > T. (3.5) 
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In the first part (a), using a standard argument and the convergence in distribution (1.7) above, the 
limit becomes as N — > oo 



max \xj\<T 3=1 

For the second part (b), the region is a union of two (not necessarily disjoint) pieces : 

(i) max \xj | = | a?x | 

3 

(ii) max |afj| = \x k \. 

3 

Note that since x\ > • • • > aifc, maxj j is either \x\\ or |xfc|. Over region (i), 



/ ni^rrfp^ lan (a<^i 



Zk<x k )< I \xi\ ai+ - +a *dVff m (£i <x u --- ,tk<x k ) 



< 



\Xi 



ja 1+ -+a fcrfp Pla„ 



(-oo,-T)U(T,oo) 

E ^(ieir i+ - +a Hx6<-T+x6>T)). 



(6 < X!) 



Similarly, 



(3.6) 

(3.7) 
(3.8) 



(3.9) 



/ n i^r^^Ki < •• < x k ) < K lan (\^r + - +ak (x^<-T + x^>t)). (s.io) 



Now from the tail estimates in Proposition L2 below, the moment (3J3) as N — > oo is equal to (3.6) 
plus a term which can be made arbitrarily small if we take T large enough. However, from Lemma 3.3, 



for T large, (3.6) is arbitrarily close to E(x" 1 • • • x c ^ k ). Thus we have proved the theorem. 



□ 



Lemma 3.3. For any k S N, and for any aj € No, 1 < j < k, 

E(|a;i| 01 ».|i fc | *) <oo, 



(3.11) 



where E is i/ie expectation with respect to the limiting distribution function F in fll.q ) awd (1.7). 
Proof. We need to show that 



/ nN a w( 



Xl, ■ ■ • ,Xfe) < OO. 



(3.12) 



Fix T > 2. We split the integral into two parts as in (3.4), ( p.q ) : (a) max.,- |xj| < T, and (b) 
max., \xj\ > T. In (a), the integral is finite. In (b), the argument yielding (3.9), (3.10) implies that 



/ TT \xj\*'dF(x u —,x k )< E(\xi\ ai+ "+*»Xx 1 >T) + E(\x k \ a -+- +a « X x k <-T). 



(3.13) 
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(Note that the additional terms corresponding to E^ lan (|£i | q1+ "' + °*X£ 1 <t) in (O) and E^ lan (|£ fe | ai +- +afc x? fc >T) 



? Plan/ 



in ( 3.10 ) are not necessary here as dF is a smooth measure.) We will prove the finiteness of the last 
two expected values for a = ai + • • • + ak- 



First, we prove that E(x"Xa;i>T) < oo for any a € No- Note that by ( |1.7| ) and fl4.8|) below, for 
xi > T , 



l-F( Xl )= lim 1 - P^ an (Ci < Xl ) = lim E^ lan ( X4>;cl ) < Ce"^ 

N^oo iV— >oo 

for some C, c > 0. In particular, we have for any a e No, 

lim xf(l - T(.ti)) = 0. 

Thus, integrating by parts, 

pOO pOO 

/ x a 1 dF{x 1 )=T a {\- F{T)) + ax^il- F{ Xl ))d Xl 
jt Jt 



3/2 



Using (1.7), Fatou's lemma and (4.8), we have 



E(x$ X xi>t) = / x a 1 dF{x 1 ) 

JT 



lim T°(l 



N 

< liminf 



vPlan^ < T)) + / J im ax a-l {1 _ pPI-n^ < 



T a (l 

Plan^a 

cT 3/2 



-^V — *oo 

pOO 

»N lan (Zi <T))+ axr\l ~ P£' Qn (£i < x 1 ))dx l 
Jt 



liminfE^ an (^ X6>T ) 

N— >oo 



< OO. 



The proof of the finiteness of the second expected value in ( p,13| ) is similar using ( |4.9| ) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



□ 



4 Tail estimates 

For the proof of Theorem [O, we need tail estimates for the (scaled) length £j. of each row, which are 



uniform in N . In this section, we obtain these tail estimates in Proposition 4.3. These estimates follow 



from the tail estimates, Proposition 4.1 for the Poissonized Plancherel measure introduced in Section 
[l], together with the de-Poissonization Lemma 4.2. 
Define 



4<*>(t) :=P f-( Afc <n) 



A e-* t 2N 



^ Nl 

N=0 



P^ lan (A fe < n). 



(4.1) 



(In j|, |j, the notation A = \fi is used. But in this paper, to avoid the confusion with the notation A 
for a partition, we use t.) The following result is proved in Section]^ using the steepest-descent method 
for RHP. Note that < 0^ fc) (t) < 1. 
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Proposition 4.1. Define x by 



2t _ x 



(4.2) 



Let k G N. There are constants C, c > and < <5o < 1 such that for large t and n, and for any fixed 
< 8 < #0j ^ e following hold true : for x > 0, 



0< l-^ fe) (t) <Ce 
0< l-^l fe) (<) <Ce- 



r 3/2 



< - < 1 - 8, 
n 

l-8< — <l, 
n 



and for x < 0, 



2t 

1 < — < 1 + 5, 

n 

, 2* 
1 + 8 < —. 

n 



0<^n\t)<Ce- c ^ a/ \ 
< t$\t) < Ce~ ct , 
We also need the following de-Poissonization lemma : 
Lemma 4.2. There exists C > such that for all sufficiently large N , 



(4.3) 
(4.4) 

(4.5) 
(4.6) 



jPlan / 

- N 



1 (A fe <n)<C0«((7V-ViV) 1 / 2 ), 
for all n € Z. 



1 - 



pPlan 

N 



(A* <n)< C(l - 4 fe) ((JV + ^) 1 / 2 )) (4.7) 



Proof. This is similar to Lemma 8.3 in || (again note that A in || satisfies A = y/i.) Indeed, the proof of 
Lemma 8.3 in || only requires the fact that < q n ,N+i < Qn,N < 1- In our case, q n ^ = P^- lan (Afc < n), 
which is clearly between and 1. The monotonicity can be found in [E4J Lemma 3.8. □ 



Now Proposition 4.1 and Lemma 4.2 imply the following uniform tail estimates. 

Proposition 4.3. Fix k € N and a € No- For a given T > 2, there are constants C, c > and No > 
such that for N > N , 



^N an (C k X (k >T) < Ce 



_ cT 3/2 



Ce" 



cJV 1 /2 



and 



<-T. 



< Ce 



_ cT 3/2 



Ce 



cN 



1/2 



(4- 



(4.9) 



Proof, (a) Bound (4.8) : Without any loss we can assume a > 0. Note that since < Aft < N, 



-2N^<i k <^^<N^. 



(4.10) 
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If T > N 5 / 6 , then the expected value in (4.8) is zero, and the bound is trivial. Thus we assume that 
T < TV 5 / 6 . Integrating by parts and using Lemma 4.2 , 



K lan (C k x ik > T ) 



S a dV p N lan (^ k < s) 



(T,JVV6) 



T a (l-F% m (fy,<T)) + 



as 



(T,JVV6) 



< CT°(1 - 03V) ((JV + ViV) 1/2 )) + C 



(r,jv5/6) 



as a - 1 (i-0W ) ((7v + Viv) 1 / 2 ))d s , 



(4.11) 



for large N, where n(s) = 2%/iV + sTV 1 / 6 . Note that since T > 2, 2(jv +^ )V2 < 1. We distinguish two 



cases 



(o o< 2(jv :ff 1/2 <i-^ 
(«) i-j< ^+ff i/2 <i 



(4.12) 
(4.13) 



where < S < 1 is a fixed constant satisfying 5 < 5q, where So appears in Proposition 4.1. 
Case (i) : For all s > T, < 2(jv+ ^ )1/2 < 2 ( N +^/ /2 <1 - S . Note that for T < s, 

v * — ' — n(s) — n \T) — ' 

/ w m . 2(jV + VAQ 1 / 2 2x^V 
n(s) > n(i ) > : ; > 



1-5 



1 - 5 



Hence from the estimate (4.3), we have 

1 _ 0« } ((iV + ^N) 1 ' 2 ) < Ce- cn ^ < Ce- cNl/2 



for T < s < TV 5 / 6 with a new constant c. Therefore, from d4.11|), we obtain 



E^ lan (4 a X« fc >T) < Ce 



cN 1 ' 2 



(4.14) 



(4.15) 



(4.16) 



Case (ii) : There is sq > T such that 



2(N+y r N) 1/2 _ 



= 1-5. We write (4.11) as 



E£ an (aX&>T) < CT a (l - ^> T) ((N + V^V) 1 / 2 )) + C / aa-Hl ~ C(i 3 (0V + ViV) 1/2 ))d S 



c 



as^l - 4& s) {{N + v^) 172 ))^ (if s > iV 5/6 , the third integral is zero) 



(T,a ) 



as a-l e -cn( S ) ds 



(4.17) 



using (Q , H , where at(s) is defined by the formula (HJ) with t = (N+^/N) 1 / 2 and n = 2^/N+sN 1 / 6 . 
As in Case (i), for s > sq, 



2(jV + VlV) 1 / 2 2VN 
n{s) > n(s ) = — - > — -, 



(4.18) 
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and hence, the last integral is less than Ce cNl/2 . For the other terms, since 



2VN ~ 2(N + VN) 1/2 = —= 4 = ; > -1 > -N 1 / 6 , (4.19) 



we have for T < s < so, 

n-2t _2y/N + sN 1 ^-2(N + VN) 1 /' 2 s-1 \ 

X{S} ~ W2)U* ~ (V7V + fiVV6)i/3 ^ ( i + _l_)i/3 " (1 + 5#73) 1/3 



as s > T > 2. Noting that sq = 2( ' N+ ^i_^ N i ( /l S)Vn ^ cqN 1 / 3 for some constant cq, we have 



,4 ' 21) 



for s < sq with some constant c > 0. Hence 



i>oo 

mPlan/va \ ^ /-.rria -cT 3/2 . r< I o— 1 — cs 3 / 2 j , /-^ -cN 1/2 

JT 



<Ce- cT3/2 +Ce" cNl/ \ 



t (4.22) 



(b) Bound : Recalling ( pL10| ), if T > 27V 1 / 3 , the expected value in (fD^) is zero and the bound 
is trivial. Thus we assume that T < 27V 1 / 3 . Integrating by parts and using Lemma 4.2, we have for 
some constants C, c > 0, 

®Tm\xt k <-T) = / (-s)<w^ an (a < s) 

J[-2N 1 / !i -T) 

= | -T|°P^ lan (£ fe < -T) + / ai-a^F^fa < s)ds 

J[-2N 1 / 3 -T) 

< <^<&-r)(W - v^V) 1/2 + C [ ai-sr^^MN - sfNf' 2 ds 

l[-2N 1 / 3 -T) 

(4.23) 



for large N, where n(s) = 2\/]V + siV 1 / 6 as before. Given T, we take iV > 1+ ^ +T!i so that for N > N , 
2 ^ N n(^T) — > 1- We distinguish two cases : 

n(-T) 

t< 2g -^ <1 + , (4 . 26) 
n(-T) 



where < 6 < 1 is a fixed constant as above. 

n(s) — n(-T) 



Case (i) : For all -2N^ - 2tfV3 < s < -T, W-f?) 1 ' 3 > 2(JV-V^ > ^ ^ h 

V / nisi — n. ( — I i — 



estimate (p6|), using T < 2N 1 / 3 , 

K lan (\C k \x ik <-T)<Ce~ c ^'\ (4.26) 
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Case (ii) : There is s > T such that 2(N ~ { ^' 2 =1 + 8. We write ( gjf ) as 



j^Plan 



(\e k \x ik <-T) < CT a ^_ T MN - VN) 1 ' 2 ) + C 



'(-so,-T) 



+ C / a(- S ) o -Vi%((iV-VAf) i/: 0efe 

J[-2JV 1 /3,- So ] 

< CT Q e- c ^^ T 'l 3/2 + C [ o(-«)°- 1 e- c l»WI 8/9 dfl 

J(-so-T) 

[~2N 1 / 3 -s a ] 



c 



(4.27) 



using (4.5), (4.6), where y(s) is defined by a; in the formula (4.2) with t = (N — \/N) 1 / 2 and n 



2*jN + sN 1 / 6 . The last integral is less than Ce cN with new constants C, c. For the other terms, 
note that since 2-/AT- 2(7V - //V) 1/2 < 2 < iV 1 / 6 and s < —T < -2, we have for -27V 1 / 3 < s < -T, 

n-2t s + \2VN -2(N - VNY^W- 1 / 6 s 

= ^ ( i + ,W/3))V — • ^ s + 1 * 2- (4 - 28) 

Thus, we obtain 



^ lan (|4 a |xa<-T) < CT a e- cT3/2 



< Ce- cT ' + Ce 



C alsj a - 1 e- c|s| ' ds + Ce 

J(-s a -T) 

3/2 . cJV l/2 



■cN 



1/2 



(4.29) 



□ 



Remark 4-4- The results ( |4.3[ ) - ( |4.6| ) /or k — 1 were given in J^J. Indeed, in ^ stronger bounds than 
(4.5), (4.6) were obtained (Lemma 7.1 (iv), (v) in Ml) : 



0<4i\t)<Ce- c W 3 , l<3<l + 5, 
0<$\t) < Ce~ ct \ 
(note A = \ft in From this, as in the Proof of Proposition 



r 2* 

1 + 6 < — 

n 



K P N laa m\x^<-T)<Ce- cT " +Ce 



-cN 



(4.30) 
(4.31) 

(4.32) 



In this paper, we only obtain the above weaker bounds (4.5) and (4.6), but they are enough for our 



purpose in proving the convergence of moments. However, we believe that the same bound (4.32) holds 
true for general k. In the next section, we indicate why we only obtain these weaker bounds (see the 



Remark before Lemma 5.1 below) 
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5 Riemann-Hilbert problem 



In this section, we will prove Proposition |4.l[ 

and (JO), note that ^ {t) > <$\t) for all k > 1 as A x > A 2 > 



For 



But for k = 1 the 



estimates (4.3), (14) were proved in |3| (Lemma 7.1 (i), (ii)), and so we have the same bounds for all 



k > 1. On the other hand, since P^ lan (A < n) < 1, we always have < <^4 fc) (t) < 1. 



The rest of this section is devoted to proving (4J3) and (4J3). We start from the formulae (see 



H |1) that (recall ^{t) = Pf ois (A fc < n)) 



det(l-S n ), 



dct(l -rS„) 



k > 1, 



with 



(5.1) 
(5.2) 

(5.3) 



This follows from, for example, Theorem 3.1 of |Uj with p + = p_ = (<, 0, 0) (see also j3(j) which states 
that for any finite subset A of Z, 



(A c {Aj — j}) - det(S(i,j)) iJeA , S(i,j) = ^for 1 ) 



i+W-j-k, 



(5.4) 



fc=i 



vPois 



denotes the Plancherel measure for Sn with N being the Poisson 



with if given by ( |5.3| ), where 

variable with mean t 2 (Poissonized Plancherel measure). Recall from (O) that (j>£\t) = Vf ms (X k < n) 



In |3j, the authors obtained the estimates (stronger than) (4J3) and (4J3) in the case k = 1 (see the 
Remark at the end of Section ||). Hence we need to prove that for any fixed k 6 N, 



^| r=1 det(l-rS n ) 
^| r=i det(l-rS n ) 



By Cauchy's theorem, we write for < e < 1, 

d k 



dr k ' r 



=1 det(l-rS n ) 



< Ce- C l-I 3/2 , 
< Ce- Ct , 

k\ 
' 27ri 



1< — < 1 + ^, 

n 



1 + 5 < 



det(l - sS„) 



2t 



ds. 



(5.5) 
(5.6) 

(5.7) 



/|s-l|=e ( S -l) fe+1 

By Remark 3 in Section 0, for ip(z) as in ( |5.3| ), S n is positive and || S n || < 1. Hence the eigenvalues 
cij of (the trace class operator) S n satisfies < cij < 1 (actually one can show that < cij < 1). For 
|s-l|=e, 



|det(l- sS„)| = JJ|1 - sa,,-| < Y[(l - (1 - e)^) = dct(l - (1 - e) S n ). 

i j 

Therefore we have 



d k 

_| r=i det(l-rS n ) 



< -det(l-(l-e)S n ). 



(5.8) 



(5.9) 
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Thus it is sufficient to prove that for fixed < r < 1 , 

det(l-rSn) < Ce~ c W a/ \ 
det(l-rSn) < Ce' ct , 



2t 

1 < — < 1 

n 

l + 5< — , 
n 



S, 



(5.10) 
(5.11) 



or by Theorem 2.1, we need to prove that for any fixed < s < 1, (note that ip given by (5.3) has no 
winding) 



(1 + *) 
(1 + 

where x is defined in Q4.2]) 



■det(l-sK n ) < Ce- c W 3/ \ 
"™det(l - sK n ) < Ce- Ct , 



2t 

1 < — < 1- 
n 

l + 6<* 
n 



(5.12) 
(5.13) 

¥ii^- (5 - 14) 

Since K n is an integrable operator, there is a naturally associated Riemann-Hilbert problem (see 
[ pl[ 0)- Let m(z; k) be the 2x2 matrix function which solves the following Riemann-Hilbert problem 
(RHP) : with contour S = {\z\ — 1}, oriented counterclockwise, 



2t 

— = 1 

n 



i(z; k) is analytic in z € C \ E, 



i + (z; k) = m_ (z; fc) 



-tiz-z- 1 ) 



sz k e t{z z 1 ) 



for z e E, 



1. 



(5.15) 



1(2; fc) — > J as z — > 00. 



Here and also in the following, the notation f+(z) (resp., /_ (z)) denotes the limiting value of lim z /^ z f(z') 
from the left (resp., right) of the contour in the direction of the orientation. In the above case, m + (z; k; s) 
(resp., m_(z;fc;s)) means \im z /^ z m(z'; k; s) with \z'\ < 1 (resp., \z'\ > 1.) In (52) of [Q, it is shown 
that 



(1 



dct(l - sK n ) = Y[ m n (0-k + 1) 



(5.16) 



where mn(0; k) is the (ll)-entry of m(z; k) evaluated at z = 0. Therefore, in order to prove (5.12) and 
(5.13), we need asymptotic results for rnn(0; k) as k, t — > 00. In the special case when s = 1, this RHP 
is algebraically equivalent to the RHP for the orthogonal polynomials on the unit circle with respect 
to the measure e^ z+z ^dz/ftiriz), whose asymptotics as k,t — > 00 was investigated in ||. The RHP 
( 5.15 ) was introduced in Q. 

There is a critical difference in the asymptotic analysis depending whether s = lor0<s<l. 
In the former case, the jump matrix in (5.15) has a upper/lower factorization, but not a lower/upper 
factorization, while in the later case, the jump matrix has both factorizations. This difference makes the 
later case much easier to analyze asymptotically. In the former case, we need a WKB type analysis which 
involves the construction of a parametrix in terms of the equilibrium measure of a certain variational 
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problem, and in the case where — > 1, the main asymptotic contribution to the RHP comes from the 
part of the circle near z = — 1. But in the later case, due to the existence of both factorization of the 
jump matrix, the RHP localizes in the limit just to two points on the circle. We refer the reader to 
for an example of the second type, and to JIt], Q for examples of the first type. 

Remark 5.1. The different analysis for s = 1 and < s < 1 gives us different estimates. Indeed, when 
s = 1, instead of ( 5.18 ) below, we have (see (6.42) of ^j) 



2t 2t 
logmn(0;fc) < + log — + 1), 



(5.17) 



which imply (4.30), (4.31). Thus, in order to obtain the better estimates (4.30), (4.31), (4.32) for general 
row £fc (see Remark at the end of Section^), we need to analysis the RHP ( 5.15| ) as s — > 1 instead of 
fixed s < 1. 



Fix < s < 1. We will prove the following estimate. 



Lemma 5.2. There are positive constants Mo and c such that when > 1 + 9l/ ^" 2/3 , we have 



logmn(0;fc) < -c\ 1 - 



(5.18) 



for large t. 



Assuming this result, we will prove (5.12) and (5.13), which completes the proof of Proposition 



4.1 



Proof of Proposition We need to prove (5.12) and (5.13). 

(a) Estimate (|5.12| ) : It is enough to prove ( pMSp for -2^ 3 n 2 / 3 5 < x < —M where M := 
Mq(1 + 25) with Mo in Lemma 5.2. As noted above, since S n is positive and || S n || < 1, its eigenvalues 
aj satisfies < aj < 1. Hence for < s < 1, by ( |5.1(j ) and Theorem 2.1 , 



< Y[ ™n(0; k + 1) = det(l - s 2 S n ) = ]J(1 - s 2 a 3 ) < 1, 



k—n 



for any n. Therefore in order to prove (5.12), it is enough to show that 

^logmn(0;fc + l) < -c\x\ 3/2 + C, 
(*) 



for some constants c, C > where the summation is over the set (note x + M < in (5.20)) 

(*) 



x + Mi -i /, 
n< k <n -T^n 1 ' 3 , 



2 i/3 



with Mi := Mo(l + S). We will show that for k in (*), 

Mq 

2 l/3 fc 2/3 



1 



< 



2t 



(5.19) 



(5.20) 



(5.21) 



(5.22) 
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Since n < k, (5.22) follows from 



2t , M 
— > 1 H — 

k ~ 2 X /3 n 2/3 ■ 



In order to show ( 5.23| ) , since 



2t _ 2t n , x x + Mi v-i 

T~^"k~^ ~ 2 1 /3„2/3K - 2 l/3 n 2/ 3 ) ' 



it is enough to checking that 

fl- 



2 1 /3 n 2/3 

which is equivalent to check that 

Mi > (l - 

But since -x < 2 1 / 3 6n 2/3 , 



>(l-Z + *)(l + . M ° 



2 l/3 n 2/3 
2 l/3 n 2/3/ JlJ0 - 



2 l/3 n 2/3 



- x + M x . , r , . Mi x Mi 



and hence (5.22) is proved. 



Now using ( p,18| ), the sum in ( 5.20 ) satisfies 



£]ogmn(0;fc + l)<£-«/l-(^) 

(*) (*) V V J 



(*) 



£+Mi 1/3 
2 l/3 



< -c 



4i 
' :| 3 
2 /2t 



3 \n 
2 (2t 



3 In 



1 - 

-1/2 
"1/2 



1 ds 

2t 



_ x+Mi „l/3\ 3/2 
2 i/3 " 



2f 

,2i 



1-^ 
2* / 



3/2- 



a; + Mi 



2^3 

3/2 / \ 3/2- 



3/2 



(5.23) 



(5.24) 



(5.25) 



(5.26) 



(5.27) 



(5.28) 



2/ 



_ n 2 / 3 (- -1 



3/2-1 



(^-X -(— ) 



where the second inequality is due to the monotonicity of the function f(y) = y/1 — Since 1 < — < 
1 + 6, we obtain ( |5.12j ). 

(b) Estimate (|5.13|) : By a similar argument as in (a), it is enough to show that 



^logmn(0;fc + 1) < -ct + C, 

(**) 



for some constants c, C > where (**) is the set 

2t 



(**) 



1 + 6 



< k < 



2t 

1 + 6/2' 



(5.29) 



(5.30) 
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For k in (**), we have 



k 1 

< -r- < 1. 



2t ~ 1 + 5/2 
Now using ( PI ), the sum in (|5.29|) satisfies 



5>g"*n(0;fc + l)<X;-c\/l- 



(**) 



(**) 



< -C4/1- 



< -ci + C. 



2f 



2/ 



1 + 5/2/ Vi + V 2 1 + < 5 



(5.31) 



(5.32) 



□ 



RHP Asymptotics and Proof of Lemma 5.2 



In the rest of this section, we prove Lemma 5.2 by asymptotic analysis of the RHP ( |5.15 ) 
Set 

k 

7?:= 2l' 



(5.33) 



Under the condition of Lemma 5.2, we have rj < 1. We denote by v(z) the jump matrix in the second 
condition of the RHP (5.15). Note that the (21)-entry of v is se 21 ^ 2 '^ where 



where log z £ R for z > 0. The critical points of this function are £ := e l9c and = £ where 



£ = -77 + iVl 



Note that — ir/2 < 9 C < ir. For z = pe l9 , consider 



F e {p) := Ref(z) = )-{p- p l ) cos (9 + 77 log p. 



Its derivative at p = 1 is 



dp 



F e (l) = cos + 77, 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



which is positive for \9\ < 9 C , and is negative for 9 C < \9\ < tt. Indeed it is easy to check that : 

(i) . When \9\ < f , F 9 (p) < for < p < 1, and F e {p) > for p > 1. 

(ii) . When | < |0| < C , F e (p) > for < p < p , F e {p) < for p < p < 1, F e (p) > for 

1 < p < p~ 1 , and Fg(p) < for p > pg" 1 . Here po is a number satisfying < po < pg, where 



pe := 



r/ — ^/ ?/ 2 — cos 2 9 



(5.38) 



and 4F 9 (p e ) = ^^(p7 1 )=0 
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Figure 1: Graph of Fg(p) for 9 = ? (case (i)), ^ (case (ii)) and 7r (case (hi)) when 77 = 15/16 

(hi). When C < \9\ < it, F g (p) > for < p < 1, and < for p > 1. 

(iv). The curve {/oe* 61 : Fg(p) — 0} crosses the circle at 90 degree. 

Typical graphs of Fg(p) is given in Figure |l| for the value r\ = 15/16 for 6 in the three different cases 
(i)-(iii). Figure || is a signature table for Re(f(z)) when r\ = 15/16. The solid curve is Re{f(z)) = 0, 
and the dotted rays represent the lines cos 8 = —r). The ± signs denote the signature of Re(f(z)) in 
each of the four components. The curve Re(f(z)) = and the lines cos# = — r\ meet on the unit circle 
at the points £ and 

2 
1.5 
1 

0.5 


-0.5 

-1 
-1.5 

-2 

-2 -1.5 -1 -0.5 0.5 1 1.5 2 

Figure 2: Curve Re(f(z)) — when r\ = 15/16 
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2 



-0.5 0.5 1 1.5 2 



Figure 3: Solid curve represents the set {pee 19 : | < 6 < 8 C } when 77 = 15/16 



Let S = Si U £ 2 , where £1 = {e l " : \6\ < 9 C }, and £ 2 = £ \ £ x . Define the function 



S(z) 



(5.39) 



which is analytic in C \ £2 : we choose the branch so that 8{z) — > 1 as z — > 00 along the positive real 
axis. Then it solves the scalar Riemann-Hilbert problem 



S + (z) = <5_(z)(l-s 2 ), ze£ 2 
d(z) — ► 1, as z —>■ 00, 



where S± has the same meaning as in the RHP ( 5.15 ). Note that 
Now set 

( 1 



mP^ (z) := m(z)S 



-a 3 



^3 



(5.40) 



(5.41) 



(5.42) 



^0 -I) 

Then (i) nrS 2 \z) is analytic in C \ £, (ii) m^(z) — > / as z — ► 00, and (iii) mV^(z) = to^t/ 2 ) (z) for 
2; G £, where 



t,( 2 )(z) = 



1-s 2 
1 



, 1 — S z ~ 



1 

1-s 2 



z G £1, 



z e £ 2 . 



(5.43) 
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Also we have 

m 11 (0)=mf 1 ) (0)(l- S 2 ) 1 -^ 

Note that the jump matrix has the factorizations 

1 -se- 2t f^8 2 {z) 
,0 1 

1 o 

-e 



1 

se 2tf(z) S -2^ 1 



s e 2tf( Z ) 6 -2 1 _ s 2 



1 -T^e- 2tf(z) 6 2 + (z) 
1 



(5.44) 



zeEi 



z g E 2 . 



(5.45) 



From (i)-(iii) in ( 5.38) ), we can take an oriented closed curve sf 3 -* surrounding and — 1, and passing 
through £ and (the solid curve in Figure ^) on which Ref(z) < except at z = i;,^" 1 . Let 




y (3) \ 
+1 \ 




a 



Figure 4: S^ 3) and Slf\ 



3 = !,••• ,6 



s (3) = s (3) u s (3) where E (3) . g the Qpen subget of E (3) satisfying | arg ( z )| < g c and s<_ 3 ) = \ £(_ 3 ). 
Similarly, we can take an oriented closed curve E^ 3 surrounding 0, but not —1 and passing through £, £ _1 
(the dashed curve in Figure H) on which Ref(z) > except at z = £, 



Again let E<_ 3) = E^ U E^ 



where is the open subset of E+ ; satisfying | arg(z)| < 8 C and E 
will be specified further below (see the third case for the estimation of \vr(z) — I\ between (5 



(3) 



.(3) _ v (3) 

+2 



\E^. The shape of E^° 



(3) 

± 
and 
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(|5.89p). Let Z ( V = SL 3 ' U Lct &>j\ 1 < 3 < 6 be °P en regions as in Figure |. Define 



i^(z) := < 



m^(z 



m^ 2 \z 

rrS 2 \z 
m^ 2 \z 

Then (i) m( 3 )(z) is analytic in C \ E^ 3 ), (ii) m^\z) — » J as 2 
z G S™, where 



i/ 3 >(z) = < 



1 







, se 



1 



1 -se- 2t ^ 2 ' 
v° 1 

,0 1 

1 ' 



^n- 2 1 



Also we have 



mn(0)=mf 1 ) (0)(l-s :i ) 



2 e ^i 3) , 







1 - se - 2t ^ 2 \ 

1 '' 


z e 












z e f4 3) 






1 1 \ 









(5.46) 



3, and mj? (z) = m^\z)v^ for 



(5.47) 



(3) 
-2- 



„2\1- 



(5.48) 



Observe that i^ 3 -* (z) — > / as t — > oo for z € E^ 3 ) \ {£, £ 1 }. Thus we expect that (z) — > I as £ — > oo. 
If this were indeed true, we would have 



log mn (0) 



sin^ 1 yr 



log(l - a 4 ) < 



2n ^ lo g(! - s 2 



1 - 



t — > oo. 



(5.49) 



But the difficulty, however, is that i/ 3 ) does not converge to / uniformly on E^ 3 ). As in [[l6), we overcome 
this difficulty by constructing a parametrix for the solution of the RHP (E^ 3 ) , v^) around the points 

Let r be a complex number satisfying < |t| < 1. Following Q, set 

f 



v:=-—\og(l-\T\'), 



a := iv. 



(5.50) 
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Define 

012 ■= TT7 v , 021 ■= 012 = • 5.51) 

rl (—a) rl (a) 

Note that 

Piafoi = v, (5.52) 

as |r(iu)| 2 = ^ginh^^) for real f ^ 0. Let D a be the parabolic-cylinder function (see, e.g. (I], ^7|) which 
solves 



^A,(C) + Q- Y + ")^<<0= t) - (••••'■•5) 



We note that -D a (C) is an entire function. 
Let the matrix 

*N=^ U !1 ! 12 ! W !V t.eCXK, (5.54) 

\^*2l(w) *22(w)/ 

be defined as follows (see p6) Section4) : for Im(w) > 0, 

*nH := e-^A^e-*"™), (5.55) 
*ia(«0 := (^ 21 )- 1 e^^£)_ (e-^ i i ( ;)-y J D_ ( e -^ < «;)Y (5.56) 

*2iM := (P^-'e-^^Dje-^ + ^D^e-^w^ (5.57) 

* 22 H := ei w D_ 4 (e-i lri a)) ) (5.58) 
and for Im(w) < 0, 

# n (w) := e** v D a (e* wi w), (5.59) 

*i 2 («0 := {M^e-^^D^^-^D^ieJ^w)^ (5.60) 

*2iH := (^j-'ei^^ACei^ + y^Ce^)), (5.61) 

*aa(«») := e^^-D-aCe***!!;), (5.62) 
The function 'J satisfies 

• ^(w) is analytic in w e C \ K. 

• For aiel. 

*+H = *_H r ~ ' T ' ~ r ^j , (5.63) 

where ^^(w) (resp., is the limit of ^(s) as s — > w with Im(s) > (resp., Im(s) < 0). 
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• As w — - > oo, 

*(w)e4* w " CT3 iy-^ CT3 = / + O^ 1 ) (5.64) 
where w~ lv denotes the branch which is analytic in C\ (—00, 0] and has modulus 1 for w S (0, 00). 

These properties can be found in |]l6f Section 4. 

Let r be the union of four rays, labeled by Tj, j = 1, • • • , 4, with the orientation as indicated in 
Figure ||. All the rays and K meet at the angle it/3. Denote the components of C \ (r U R) by Qj, 




Figure 5: Tj and £!j 

j = 1, • • • , 6 as in Figure @. Define H(w), analytic in C \ T, by 





ff N - ( "'I' )(io )*We'" 2 ^- ( ? ( f ° - x ) , 



where a(£) is 



a(0 - e 



-^-1)^(1 - 772)1/4 



(5.65) 



(5.66) 
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and 4>(w) is defined by 



<j>(w) 



1 —re 2™ ui 
,0 1 
1 

-re~i iw2 w 2il/ I 



_r „— hiw 1 liv i 

_ T p 6 % W I . 

i?w 2 „,,— 2if ' 



1 T^FF e '™ w 
ii 1 



w G fii, 



w 6 fig 



to G fi 3 , 



ui g r24, 



to e si 2 , ^5- 



(5.67) 



Then by recalling that to W<T3 is analytic in C \ (— oo,0], one can directly check that H + (w) 
H-(w)vh(w) for w G r, where Vh(w) is given by 



1 -Ta(£)~ 2 w~ 2w e2™ 



1 



:= < 



ja{0 w z 
1 



2„„2ii< a -^iw 2 



1. 



Also, from (5.64), we have 



.T=pFpa(0 w e 
1 



H(w) = I + 0(w 1 ), asw-^oo. 



w g r x , 



u; g r 4 , 



w g r 2 , 



ui g r 3 . 



(5.68) 



(5.69) 



As |a(£)| = e _I/0c , — 7r/2 < 8 C < n, we see that the error term 0(w~ 1 ) in ( 5.69 ) is uniform for ^ > 1. 
Similarly, |iJ(w)| is uniformly bounded in the w plane for ^ > 1. 
Define the map 



z ^ w(z) := V2t(l - rff^iC 1 ^ - £). 



(5.70) 



It maps £ to 0, and the tangent line Lj to the unit circle £ at £, to the real line as in Figure §. Let Oe, 
Oj be the disjoint sets {z : \z — £| < g}, {z : \z — £| < g}, respectively, where g is defined by 



1 i Mo < 2t < 1 , : 

1 T 91/31.2/3 ^ h ^ 1 ' "> 



l + 5< f • 



(5.71) 
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w(z) 



Figure 6: map z i— » w(z) 



The (small) parameter < e < 1 will be specified below (see ( 5.105| ) below). We note that one may 
choose the curves in E' 3 ) above so that in 0£, Oj, they are straight lines which map under z i— > w{z) 
to (finite subsets of ) the rays Tj, j = 1,- • • ,4, E^ flO^ T 4 , E^ flO^ T 2 , E^ flO e ^ Ti, 

ESSn. 



— * Ts, and similarly for the neighborhood of Oj. For r — s, we define 



H(w(z)), 
H(w(z)), 



(5.72) 



zeC\(O c uO ? ). 



Let E B := E< 3 ) U<90 ? U<90t as in Figure H where and Or are oriented counterclockwise. Clearly : 




Figure 7: E fl := E^ 3 ) U d£> 5 U d£>r 



solves a RHP on Ej? : m p (z) is analytic in C \ E^ 3 \ m p (z) — > / as z — > oo, and m p+ (z) = m p -(z)v p (z) 
for z g Efl for a suitable jump matrix v p . Set -R(z) := m^ 3 ' (z)m p (z)~ 1 . Then R+(z) = R-{z)vr{z) for 
zeEjj where = mp-u^ 3 ^" 1 /?!"^. Now we estimate \vr(z) — I\. 
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For z € dO ( , we have from fls/Fl] ), when 1 + 21/ f° 2/3 < f < 1 + 5, 



\w(z)\=e(l-? 1 2 f^V2i 



> e 



3/4 



(5.73) 



> 



eM, 



3/4 



2V4(1 + <J)l/4- 



When 1 + 6 < f , 



;(z)| = v^(l-7? 2 )e> f^J eV2l. 



1/4 



Thus if we have taken Mq large, and t is large, we have for z 6 <90£, from ( 5.69 ), 

1 



m p+ (z) = I + O 



iin(M 3/4 ,Vt) 



But as = I on <90£, vr(z) = v p (z) 1 = m p+ (z) 1 , and hence 

C 



\\vr(z) -I\\ L °°(dOs) 



< 



min(M 3/4 , Vt) 



z e dO € . 



(5.74) 



(5.75) 



(5.76) 



We are using here the standard fact that if deti> p = 1, then detmp = 1. Similarly, we have the 
same estimate ( 5.76 ) on v R (z) for z g dOj. 

For z e s (3) no f , since m p and m p 1 are uniformly bounded, \v r(z) — 1\ < C\v^ (z)v p (z) 1 — I\. 
For z e £^ n O s , by (|5~4% , ( ^68|) , 



^(zKW- 1 -/! < |4 3) (z)-(^) 21 Mz))| 



5|e 2t/(z) (5" 2 (z) - a(^) 2 u)(^) 2 ^e-i iw W 2 1 =: s|A|. 



(5.77) 



Setting u := i£ (0 — £), we have 



£(l-iu)-£- x 

U-rv 

(e 2t,l(u) j(w) - l)e 2 * ( ^-«~ 1 )+^°s«)- ti ( 1 -'? 2 ) 1/2tl2 



(5.78) 



where 



1 . 1 



) + r) log(l - iu) + - ?7 2 ) 1/2 7 



(5.79) 



(-g^ + ^vT^V + cV), 
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and 



JM = I J^fsf = 1 + 0( ^» " 1 + 0< 7^? ) ' ,5 - 80) 



Also, as I , " I < ce for z £ Oc, we have 



< c|u| 3 , (5.81) 
= l + 0(_==) = l + 0(e). (5.82) 



On the other hand, for z £ E^j n O^, 



Rei-ityjl-jfu 2 ) < -cty/l-r} 2 \u\ 2 , c = cos ~ > 0. (5.83) 

6 



Therefore, we obtain 



|A| < C\{e 2th(u) - 1) + e 2th{u) (j{u) - l)\ e - ct V^\ u \ 2 

< C(\2th(u)\ + \j(u) - l^ e -ctV^f\u\ 2 +\Re{2th(u))\ 



< 



C (t\u\ 3 + e)e- ct V^f\-\ 2 (5 - 84) 
C 

< —= + Ce, 

" V2t{l - r? 2 ) 3 / 4 

where for the last inequality, we have used the fact that \x 3 e~ x | is uniformly bounded for i£l. 
Now for 1 + 2l/ Mo, /3 < f < 1 + 5, we have 

V2i(l-^>VS(l (|)" 4 (^(f -l,) 3 ' 4 ,^^) 3 ' 4 , ,5.85, 



and for 1 + 5 < ^, we have 



>/2t(l - r? 2 ) 3 / 4 > ^(^) 



3/4 



(5.86) 



Thus we obtain (recall (5.77)) 



\\vR-I\\r^rS°}noA < Cs W&\\r,~<rS°}nnA < . ,.Sm ^ + ^ ( 5 - 87 ) 



which is small if we take Mo,t large and e small. For other parts of E' 3 ) nO^, by a similar 
argument, we obtain the same estimate. By the symmetry m^(z) — m( 3 )(z) and m v {z) = m p (z), 
we obtain the same estimate for E^ 3 ) n Oj. 

Let O := C ? U (9^. For z E E^ n (C \ 0), «k( z ) = f (3) (^)- Thus we need an estimate for 
1% (#) = se 2t -^'<5 _2 (z). Since |<5(z)l = e _lye where 9 = arg(^ z =), \S(z)\ and |<5 _1 (2:)| are uniformly 
bounded. When 1 + 5 < 4^, disf(E_] D 0°, is uniformly bounded below. From this fact, 
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one can check that we can take £_j so that Re(f(zj) < — Co(e) for z S S_jflO for some constant 
Co(e) > depending on e. Hence we have 

2+ 

IK - /|l L . (E wr*n < Ce ~ co(e)t > 1 + «^ p ( 5 - 88 ) 

On the other hand, when 1 + 2l/ ^ ] 2/3 < f < 1 + (5, we take £^ = {/5(6»)e ie : |0| < 6< c } such that 

(i) For z £ S^f] with ^ < |ar<7(z)| < 6 C , is a pair of straight lines which meet the unit 
circle at £ and £, respectively, with angle 7r/3. 

(ii) For pe i9 E SL 3 J with |ar 5 (z)| < ?f, Re{f{pe ie )) < Re(f{p'e ie )) for p < p' < 1. Also is 
an increasing function for < 6 < and is a decreasing function for — ^ < 8 < 0. 

(Here the precise value ^ is of no importance : any angle between ir/2 and it will do.) Condition 
(ii) can be achieved by choosing S_| always to be above the curve {pee 19 ■ § < < 9 C } (recall 
( 5.38 ) and Figures 0, |||). Condition (i) can be achieved as the curve {pee 19 ■ § < < &c} 
crosses the unit circle at 90 degree (see Figure |3|) For z in (i) satisfying arg(z) > 0, we have 
z = £(1 — ire~% 1 ) for some real r > 0. We note that r < < For such z, we have (recall 

( B ) 

i?e(/(z)) = A(r)^l- V 2 + B(r)r], (5.89) 

where 

A(r) = , fi(r) = ^(r-VTW-VZr) 1 _ 

4(l-V3r + r 2 ) W 4(l~V3r + r 2 ) 2 8V ; V 7 

One can easily check that A(?') < for < r < v3 and B{r) < for < r < Thus for z in 
(i) satisfying arg(z) > 0, we have for some c > 0, 



Re(f{z)) < A(r)y/l -if < -cr 2 ^l-7f. (5.91) 

For z in (ii), note first that for fixed < p < 1, Re(f(pe l6 ) is an increasing function in < 9 < it. 
Let Zf, be the point on sf 3 ^ satisfying arg(z) = 4p Thus together with the condition (ii), we 
obtain for z in (ii) satisfying arg(z) > 0, 



Re(f(z)) < Re(f(z b )) < -c\z b - t\'Wl-V 2 < ~c\z - £| 2 \A-^- (5.92) 

i 1 fo\ c 

Here the second inequality follows from (5.91). Thus we have for z € T,_{ (1 with arg(z) > 0, 



\v R {z)-I\ < { , " (5.93) 

Ce--l^lVW, l + 5I7 ^<f <! + & 

By symmetry, we have similar estimates for z £ with arg(z) < 0. Since |z — £| > g for z G O , 
the above estimates imply in particular that 

7-11 ^ >/ ' 2l/3fe2/3 ^ k J- -T 




3G 



For £^ n (C \ 0), by the symmetry Re(f(pe l9 )) = Re{f{p- 1 e l0 ), we have the 



same estimate. 



(3) ("3") — 

Also by a similar argument, we obtain a similar estimate for (£_2 U £+2) ^ \ C)- 
As usual, define an operator on L 2 (T,r), 

C VR (f) = C-(f(vR-I)) (5.95) 

where C_ is the Cauchy operator 

/(*) 



(C_/)(z) = lim 



2™ 



— ds, z 6 z' on the — side of £_r. 



(5.96) 



As the Cauchy operator is scale invariant, C_ is bounded from L 2 (Y,r) — > i 2 (Sfl) uniformly for 

f > 1 + 2 i/3fc2/3 ■ and we have II^hII < I for *> M ° sufficiently large by ( j576| ), j^87|) and 
Hence 1 — C!^ is invertible. By standard facts in Riemann-Hilbert theory (see |ll], ^|), the solution 
R(z) to the RHP (Y,r,vr) is given by 

f (I+(l- C VR )- l C VR I)(v R -I)(s) 



R(z) 



1 

2~7ri 



s — z 



(5.97) 



As TOp(O) = J, we have m$(Q) = i?n(0). By using dist(0,S fi ) > 0, ||(1 - )^ 1 1| < c, and ||C_|| < c, 
we have 

l™( 3 )/n\ 



(5.98) 



as — i]|i«> is bounded. We estimate \\vr — in each part of S fl . First, for dO and E^nO, 
since the length of the contour is of order g, we obtain by (5.76), (5.87) 

\\vn-I\\ LH ^<C e [ ' 



< 


c 


VR 


-I\ 




r) +c||(l - C VR ) 1 C Vr I\\ L 2 {Tjr) \ 


Ufl--f||L=»(BH) 


< 


c 


VR 


-I\ 


V- + 


cHi-c^y'h^WC-ivR- 


I)\\l4vr-I\\l* 


< 


c 


Vr 


-I\ 


IA + 


c\\ v R- T \\h 




< 


c 


Vr 


~I\ 


IA + 


c\\vr-I\\ l ~\\vr-I\\ L i 




< 


c 


VR 


~I\ 


L!(S 


b) 





(5.99) 



When 1 + 



2 l/3 fe 2/3 



< f < l + (5,by(^7ll), g^e x /T~^2. When 1 + ^ < f, V /T _ ^2 = yi - WMF > 
C, and hence we have g — e < cy 1 — r/ 2 . Thus in both cases, we obtain 

W-I\\ L ^o)<CV^{ L (5.100) 

Now we compute - ^|| L i (SRn o c )- We nrst focus on S -i n ^ n {im(z) > 0} When 1 + 5 < by 



(3.93), 



C 



(5.101) 
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for large t as ^/l - rf > C in this case. When 1 + 21/ fS 2/3 < t < 1 + S > from (|5.93|) 



K / llLi(s ( _ 3 ;nO c n{/m( 2 )>o}) - / (3) 



SLjflO n{Jra(a)>0} 
. oo 

< C / e^V 1 ^ 2 ^^ 
C 



Ce -c tx /rV| Z - 5 | 2 | dz | 



< 



,-cty/l-r,* 



But since, for 1 + 2l/ f° 2/3 < f < 1 + S, 



1/6 



CM 



1/3 



we obtain 



C 



IK - J ILi(E(«n-ern{/mW>o}) ^ TT374 e V 7 !^/ 



(5.102) 



(5.103) 



(5.104) 



By a similar computation, we obtain the same estimate for the other parts of d O . Thus if we take 
e small, and then take M , t large, we obtain by ( |5.98 ), (5.99) and ( 5.104] ), 



(5.105) 



with a constant a > which can be taken to be arbitrarily small. Therefore, from (5.48), ( 5. 98] ), using 
( 5.105 ), we obtain (note ( 5.49 )) for large t, 

logrnn(O) = logm^O) + (1 - — )log(l - s 2 ) < a y/l - if - cy/l - if < -C^/l-if, (5.106) 



for some C > 0, which is (5.18). 



6 Multi-Painleve Functions 



In this section we will show that the multi-interval case considered in Theorem 2.12 is related to new 
classes of "multi-Painleve function" . As we will see, these functions describe the interaction of solutions 
of Painleve equations in a way which is strongly reminiscent of the interaction of classical solitons. 
We suggest the name "Painlevetons" or simply "P-tons" for these functions. In this section we only 
illustrate a few of the properties of P-tons. The general theory will be developed in a subsequent paper 
together with Alexander Its. 



From Theorem 2.12, in the k interval case, 



(sj - sj-i) K nj (z, w) = — s-y , (6.1) 



w 
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where 



/ = (/o, • • • , fkf = {s k , (si - s )ip{z)z ni , • ■-,(«*- s/c-iV(z)z™ fc ) T , 
.9 = (.90, ■ • • ,.9*) T - -{<p{z)z n -)-\ -[<p{z)z n «y l ) T . 



(6.2) 
(6.3) 



Thus by the integrable operator theory |2l], the associated jump matrix vonE = {\z\ = 1} has the 
form 



v = I — 2nifg = 



( l-s k Sfe(^ ni ) _1 
-(si - soVz" 1 



l<p,g<A; 



(6.4) 



For purposes of illustration, we will only consider the case when k = 2, 

/ l-s 2 s^z"!)- 1 s 2 (^)-^ 

^-(S 2 -Sl)^" 2 (S 2 - Sl)^™ 2 -™ 1 1 + S 2 - Sl y 



(' = V 



(3) 



(6.5) 



and (/? = e l ( z z ) as in Introduction. Observe now that when s% = 0, the jump matrix takes the form 

/ 1-82 S 2 (<fiZ n i)~ 1 S 2 (^)" lN \ 

1 



v = v 



(3) 



Let m (3) be the solution of the 3 x 3 RHP 



s 2 z- 



(6.6) 



i + s 2 y 



mv' = m s _'v y "> , z£S, 
m^ 3 ) — » J as z — > oo. 



(6.7) 



But it is clear that the 2x2 matrix mS 2 ^ constructed from as follows, 



777. 



/ (3) (3)^ 

(3) (3) 
« TO 33, 



(6. 



solves the RHP 



(2) (2) / 1 - S 2 S 2 (^Z rl2 ) ] 

^-s 2 ^ n2 1 + s 2 
—* I asz-^ oo, 



(6.9) 
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which is an RHP which is algebraically equivalent to the RHP for Painlcvc III (PHI) which occurred 
in 01 : set 



,( 2 ) = 



,( 2 ) = 





V 1 + S2 / 

VT+~S2 

1 



,(2) 







,(2) \/T+S2 







\z\ < 1, 



Ul > l. 



o 



vT+la, 



Then m^ 2 ) solves the RHP 



.(2) _(2) / S 2 (^" 2 )- 1 

= m_ 

-S 2 (/3Z™ 2 1 



as z — > oo, 



which is the RHP for PHI considered in 18 . On the other hand, if s\ = s 2 = s, then 

/ 1-s s^z" 1 )- 1 s(^ n2 ) _1 \ 



,0) = 



Now 



-S(fZ 





,(2) 



1 + S 





' (3) (3P 

™il "H 2 

(3) (3) 

m 21 m 22 



/ 



solves the RHP 



™(2) (2) 



1-s s(ipz n2 Y 1 
m(2) — > J as z — * oo, 



which again is the (equivalent) RHP for PHI. Also if we set n\ = ri 2 = n, 



,(3) = 



I l-*a M^T 1 s 2 (^™)- 1 \ 

— Sl^SZ™ 1 + Si Si 

y-(s 2 -si)^z" (s 2 - si) l + s 2 -Siy 



Conjugating the solution of the RHP associated with by 










0^ 









o\ 







1 










1 







1° 


1 


V 




v° 


1 


V 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



(6.15) 



(6.16) 
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wv 



find that -> I as z 



It follows that necessarily 



and hence 



oo, and solves a RHP with jump matrix 

' l-s 2 s 2 (<pz n )- 1 ^ 
v {3) = - Sl ipz n 1 Si 

^-s 2 tpz n 1 + s 2 ) 

(mg mg) T = (0 1 Of 



,(2) 



;(3) 



"11 "'13 

•(3) ~(3) 



'33 



solves the RHP 



(2) (2) / 1 - S 2 52(^2") 1 

\-S 2 lfiZ n 1 + s 2 
— > I as z — * oo, 



which is again the (equivalent) RHP for PHI. 

The analogy with solitons is particularly clear if we consider in the edge scaling limit, 



2r + ^r 1 / 3 , 



j = 1,2 ; h<t 2 



as t — * oo. Then 



,(3)/ 



2iu , 



on the real line, where 



I 1 - s 2 
\-(s 2 - s\)e~ 



-u +tjU, 



„2i8. 



s 2 e" 
^2i(Si-e 2 ) 



s 2 e 

1 + si sie 
{82 - Sl )e 2i ^- ^ l + s 2 - Sl J 



J = 1,2. 



(6.17) 



(6.18) 



(6.19) 



(6.20) 



.21) 



(6.22) 



(6.23) 



In addition to varying si, S2, we can now vary ti, £2- In particular, we can follow the trajectory of the 
solution of the RHP as t 2 moves from t\ to oo. As t 2 — > t%, the solution becomes Painleve II (PII) and 
as t 2 — > oo, it gives to another solution of PII, but now with a phase shift (see |Q). It is this behavior 
of P-tons, in particular, that is reminiscent of soliton interactions. 



7 Colored permutations 

First, the definition: 

Let 7r be an m-colored permutation (see, e.g., and assume the colors are indexed by 0, 1, ... m — 

1. Let 5 be a subsequence of length I of 7r which is a union of monochromatic increasing subsequences; 
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let ki be the number of these sequences having color i, and set k — fej. Note that the monochromatic 
increasing subsequences may be empty, but the color of empty subsequences still matters. We assign 
to S the following score: 

na+ ( k V) + £ (^-4 hi+ 2 1 )) (7.i) 

^ ' 0<i<m-l ^ ' 

Now, let Ik (t) be the maximum score over all unions of k monochromatic increasing subsequences (note 
1 (tt) — 0). We then define 

A fe (7r) := Z fc (7r) - h-ifr). (7.2) 

Lemma 7.1. Let Xj? (it) be the partition associated to just the i-colored subsequence ofn. Then Afc(7r) — 
k is simply the kth largest of the numbers m(\j (it) — j) + i- Moreover, if n has length n, then Afc(7r) 
is a partition of ran. 



Proof. Fix a composition ki, and consider the largest score associated to that composition. Clearly, we 
can maximize the score for each color independently; we thus obtain: 




Now, for a fixed value of k, this is clearly maximized when the values m(\j (n) — j) + i occurring in 
the sum are chosen to be as large as possible. Plugging the resulting value of ife(7r) into the formula for 
Afc(7r), we obtain the first claim. 

Note that the numbers m(A~- (n) — j) + i are all different (the congruence class modulo m depends 
on the color, and the numbers are distinct within a given color) . Furthermore, we readily verify that for 
each congruence class, the number of negative numbers not occurring in the set is equal to the number 
of nonnegative numbers occurring in the set. We thus conclude that Afc(7r) is indeed a partition. It 
remains to verify that ^2 k Xk(^) — run; in other words, lk(^) = mn for k sufficiently large. Choose k 
such that 7r is a union of k increasing subsequences, and consider l m k{^)- We readily verify that the 
term 

E (iki-m( ki + 1 )) (7.4) 

0<i<m-l ^ ' 

is maximized when all ki are equal to k, and thus the optimal score differs from mn by 

( m \ +1 )+ E Hk-m( k+ 2 1 ))=0. (7.5) 

^ " ' 0<i<m~l ^ ' 

□ 

Remark 7.2. An alternate approach is to define Xk(^) via the Schensted correspondence for rim-hook 
permutations given in fi34j , at which point the lemma follows immediately. The fact that the rim-hook 
correspondence splits into m ordinary correspondences gives the increasing subsequence interpretation 
above. 
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Now, suppose we choose n randomly according to a Poisson law of mean mt 2 , and then choose an 
m-colored permutation of length n at random. Equivalently, take m independent Poisson processes in 
the unit square (one for each color), and convert the resulting point set to a colored permutation. We 
thus see that the resulting random partitions are independent, and are all distributed according 

to the law for ordinary permutations. In particular, we obtain the following correlation kernel: 



S (m) (a,6) = ^2(<f 1 )(a+k)/ m l P(b+k)/ m (7.6) 

k>l 

where 

<p(z) = e^*-*~^. (7.7) 



(Recall from Corollary 2.1C that ip a and ((p^ 1 ) a are for a non- integral.) 

Now by using Corollary 2.10 and Theorem [O] for the convergence of moments for the ordinary 
permutations, we obtain the convergence of moments for A^'s in the colored permutation setting. More 
precisely, as in (3.1), there is a limiting distribution F color ( m ) such that 



ii^ m 2/3( mi V)l/6 



where ]g™' or ( m ) denotes the expectation with respect to the natural counting measure on the col- 
ored permutations (see |||), and E color(m) is the expectation with respect to F color( - m \ The function 
F color ( m > (xi, ■ ■ ■ ,Xk) has the following meaning in terms of CUE. Take m random GUE matrices of 
size N at random, then superimpose their eigenvalues. We denote the largest of those superimposed 
numbers by Zi(iV), the second largest by z 2 (iV), and so on. Then F color{ - m \x 1 , ■■■ ,x k ) is the limiting 
distribution of Z\, • • ■ , z& as N — > oo, after appropriate centering and scaling. 

A number of other statistical systems which are currently of interest can also be analyzed by the 
methods of this paper. In particular, we have in mind the random word problem ]35|, [24|, ^2, 23 , certain 
2-dimensional growth models J25| , and also the so-called "digital boiling model" |2(J . 

For example, in the growth model considered by Johansson in let a = Uj =1 (Tj be a union of 
k disjoint increasing paths <jj in the model. Let LW(cr) be the sum of the lengths of the paths uj, 
and let = max CT L^(a). We define = — L^ k x \ The joint probability distribution for 
Ai, • • • , Aft can be obtained |25| by various differentiations of det(l + s jX[n j .n :j -i) S) w hh respect 

to si, ■ ■ ■ , Sk as in (1.12) with <p now given by (p(z) — (1 + ^/qz) M (1 + y /qz~ 1 )~ N . But now by Theorem 
2.12, det(l + y^ 8jX[n,j .n.)~i) S) can be expressed in terms of the determinant of an integrable operator 
as in (2.54). This opens up the possibility for the asymptotic analysis of the convergence of moments 
for the joint distribution. However, the associated RHP has a new feature, namely the weight function 
is non-real, which has not yet been addressed in general (however, see p7| ). There are similar formulae 
for random words and digital boiling. 
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